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SUMMARY 

Hypersonic flows over cones and straight biconic configurations are calculated for 
a wide range of free stream conditions in which the gas behind the shock is treated as 
perfect. Effect of angle of attack and nose bluntness on these slender bodies in air is 
investigated extensively. The numerical procedures are based on the solution of complete 
Navier-Stokes equations at the nose section and parabolized Navier Stokes equations 
further downstream. The flow field variables and surface quantities show significant 
differences when the angle of attack and nose bluntness are varied. The complete flow 
field is thoroughly analyzed with respect to velocity, temperature, pressure and entropy 
profiles. The post-shock flow field is studied in detail from the contour plots of Mach 
number, density, pressure and temperature. Flow separation is observed on the leeward 

0 O 

plane for an on-axis, 12.84 / 7 (fore-cone and aft-cone angle) biconic geometry at 

0 

12 angle of attack. Also, the windward and leeward heating rates for the fore-cone 
section decreases by a factor of four and by a factor of five respectively, when the nose 
bluntness is increased by an order of magnitude. The effect of nose bluntness for slender 
cone persists as far as 200 nose radii downstream. 
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Chapter 1 
INTRODUCTION 


The field of hypersonic aerodynamics has been studied extensively in recent years 
with a resurgence of interest. The main reason for this motivation is the currently 
envisaged National Aerospace plane (NASP) which will cruise at hypersonic speeds. The 
propulsion system of such hypersonic flight vehicle is dominated by complex interactions 
of multiple shocks with the viscous flow. A better understanding of the complex flow 
field in different regions of the aerospace vehicle and determination of performance 
characteristics is essential for vehicle design. The renewed interest in hypersonic flows 
is also due to possibility of manned exploration of the planets. Mars, because of its close 
vicinity to earth and its physical characteristics makes it the first choice for such a mission. 
After adequate knowledge about Mars atmosphere being available, the aerocapture and 
aerobraking have become promising techniques in planetary missions and Earth orbital 
transfer applications. With the use of these techniques there is a significant increase in the 
payload weight and reduced transit times as compared to a pure rocket propulsion system. 

The aerobraking [1]* technique uses the drag during successive passes through the 
upper atmosphere of the planet to circularize a highly elliptical orbit. The aerocapture 
represents a transfer of the vehicle into a closed stable orbit from a fly-by trajectory. 
The primary design factors for these transatmospheric vehicles are high lift to drag ratio, 
low hypersonic ballistic coefficient, and high volumetric efficiency. Straight and bent 
biconics are being considered as the relevant configurations for these vehicles. Also, 
biconic shapes would fit nicely into the shuttle orbiter bay. 


The numbers in brackets indicate references. 
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To avoid a conservative design involving larger and heavier control characteristics, 
thus reducing scientific return, an accurate database on vehicle performance during entry 
is needed. Because wind tunnels are not capable of duplicating the full flight trajectory 
of a planetary entry vehicle, numerical techniques must be relied upon heavily to furnish 
this database, that is, to furnish the flow environment surrounding the vehicle during 
entry. With recent advancements in Computational Fluid Dynamics (CFD), a wide 
variety of complicated flow fields in different regions of the transatmospheric vehicle 
can be analyzed. Numerical modelling of the flow field has proven to be a valuable tool 
for getting better insight into the complex nature of these flows. 

For many hypersonic applications, forebody configurations are relatively simple 
and can be modelled by long slender bodies of revolution with a blunted nose. In the 
analysis of the vehicle flow field , the forebody flow provides the initial conditions. It 
becomes inevitable to accurately predict the forebody flow field by incorporating realistic 
flow conditions and geometrical modelling. There are two important features which have 
significant effects on the entire flow field, one is the angle of attack (AOA) and the 
second is the nose bluntness. 

One of the important fields of investigation in hypersonic aerodynamics is the 
determination of heat transfer to the body surface as a function of body shape and flight 
conditions. This problem is of practical interest because the aerodynamic heating at 
hypersonic speeds is very large and represents the major limiting factor for some of 
the useful flight conditions. Usually, all hypersonic vehicles have blunt noses or blunt 
leading edges, in order to reduce (to acceptable values) the heat transfer in the stagnation 
regions. Blunt leading edges, however, result in a drastic increase in drag. With the 
advancement in material technology, it is possible to reduce the leading edge thickness 
without encountering intolerable heating rates in the vicinity of the nose. 

The presence of blunt noses or leading edges and, therefore, of detached and highly 
curved shocks introduces the necessity of considering in the heat transfer investigations 
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new parameters related to the entropy gradients existing in the flow. The presence of 
large entropy gradients in the flow field can affect the heat transfer on the body and the 
transition of boundary layers from laminar to turbulent. The extent to which this bluntness 
influences the development of flow downstream is not known precisely for hypersonic 
flows. In general, such a calculation is very difficult because of the various intervening 
effects such as pressure interaction, vorticity interaction, and curvature parameters. 

When the free stream is at an angle of incidence, the wall quantities and flow 
field variables are affected significantly. The circumferential shock shape gets displaced 
towards the leeward side and no longer remains concentric with the conical body. The 
surface pressure, wall heat transfer and skin friction show higher values on the windward 
meridian as compared to the leeward meridian. The circumferential distribution of surface 
quantities show maximum values on the most windward rays and minimum on the most 
leeward rays. At AO A exceeding the fore cone half angle, there exists an imbedded 
shock within the shock layer on the leeward side. At high AOA, the flow on the leeward 
side separates in a circumferential mode, forming two symmetrical, supersonic, counter- 
rotating longitudinal vortices. 

To calculate the flow over blunted slender bodies, two sets of equations are utilized 
in the present investigation. For the nose section, the use of full Navier-Stokes (NS) 
equations is necessary due to the presence of a relatively large subsonic region, but 
immediately downstream of the nose region on the body, the flow becomes predominantly 
supersonic and the space marching parabolized Navier-Stokes equations (PNS) are used. 

Recently, Singh and Tlwari [2] have done an extensive literature survey on the 

nose bluntness, AOA and viscous interaction effects on slender bodies at hypersonic 

speeds. However a few important numerical and experimental works related to the present 

0 

study are reviewed here. Tracy [3] conducted an experimental study of a 10 semi- apex 

O 

circular cone in air at Mach 8 and at yaw angles up to 24 . The basic measurements 
consisted of surface pressure, heat transfer and a complete survey of the Pitot pressure 
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in the flow field. The probe data delineated the boundary between viscous and inviscid 
flow and determined the shape of the outer shock as well as imbedded shocks which 
appeared in the flow field at large yaw. Lubard and Helliwell [4] investigated the flow 
field on cones at high angles of attack for hypersonic laminar case. They developed an 
iterative implicit finite difference technique to solve the three-dimensional PNS. The new 
theoretical approach was validated by comparing the numerical results with experimental 

O O 

data obtained by Tracy [3] on a sharp 10 half angle cone at a=12 with Mach number 
8 and by Stetson [5] on a sharp 5.6° half angle cone at a=8° with Mach number of 
14.2. A good agreement was found between experimental and numerical work. Stetson 
[6] studied the effect of bluntness and AOA on boundary layer transition on cones and 
biconic configurations. The rearward displacement of transition due to bluntness was 
found to be quite sensitive to free stream Mach number and entropy swallowing distance. 
Transition location was also found to be sensitive to small changes in AOA. 

Holden [7] conducted an experimental study to examine the separate and combined 
effects of surface roughness, entropy swallowing, and boundary layer transition on 
the distribution of heat transfer and skin friction over a series of biconic nose tip 
configurations and on sharp and blunted slender cones in high speed flow over highly 
cooled surfaces. In this investigation, emphasis was placed on examining the effects 
of free stream Mach number, Reynolds number, and wall to freestream stagnation 
temperature ratio, in addition to exploring those introduced by roughness height and 
model configuration. Miller and Gnoffo [8] conducted an experimental investigation 
of hypersonic flows over biconics at incidence. Pressure distributions, aerodynamic 
coefficients, shock shapes, and oil flow patterns were measured on spherically blunted, 

0 0 o 

12.87 / 7 on- axis and bent biconics (fore cone bent 7 upwards relative to aft-cone) 
at Mach 6 in air for various AOA’s. The experimental results were compared with the 
inviscid Stein code and PNS code of Vigneron et al. [9]. Miller et al. [10] did an 
experimental study to determine the laminar heat transfer distributions on spherically 
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blunted on-axis and bent biconics at hypersonic-hypervelocity flow conditions in the 
atmospheres of helium, hydrogen, air and carbon dioxide at different AOA’s. The effect 
of nose bend and real gas flow phenomena on heating distributions for biconics was also 
studied. Experimental results were compared with a three-dimensional PNS code. Gnoffo 
[11] modified the PNS code of Vigneron et al. [9] for straight and bent biconic geometries. 
Real gas effects were included by coupling the vectorized chemical equilibrium PNS 
code of Kumar et al. [12] with a variable-effective-gamma capability. With the variable- 
effective-gamma option, gamma has been defined as the ratio of enthalpy to internal 
energy of the gas mixture referenced to 0 K, to replace the usual definition as the ratio 
of specific heats. Some recent development in this field can be found in [13-18]. 

Moskovitz et al. [13, 18] conducted an experimental investigation of the effects 
of surface perturbations on the asymmetric flow past a slender body for laminar flow 
conditions. The perturbations included variations in model tip sharpness and roughness 
as well as discrete surface perturbations represented by cylindrical beads of varying 

O 

height and diameter attached to a cone/cylinder model having 10 semi-apex angle and 
3.0 caliber tangent ogive model both having base diameters of 3.5 inches. The tests were 

O 

conducted at 40 AOA. They concluded that local surface perturbations could reverse the 
direction of flow asymmetry of the basic model if they were large enough relative to the 
body radius where they occurred. The perturbations were most effective when located at 

. O 0 

an azimuthal position between 120 and 160 from the windward ray. 

Prabhu and Tannehill [14] developed a new PNS code to compute the hypersonic 
laminar flow of a multicomponent, chemically reacting mixture of thermally perfect gases 
over two dimensional and axisymmetric bodies. The new PNS code solves the gas 
dynamic and species conservation equations in a coupled manner using a non-iterative, 
implicit, space marching, finite difference method. Ericsson [15] studied the effect of 
nose bluntness and cone angle on boundary layer transition for slender cone geometry. 
He simulated the combined effect of these two parameters on transition analytically. 
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He also showed how the transition dependence on the nose bluntness induced inviscid 
shear layer changes when moving from the entropy swallowing (small bluntness) to the 
entropy impingement (large bluntness) region. Malik et al. [16] studied the effect of 

O 

nose bluntness on boundary layer instability for Mach 8 flow past a 7 semi-vertex cone. 
In order to account for the variable entropy effects associated with the bow shock, the 
basic flow was computed using PNS equations. Linear stability analysis of the basic 
flow revealed that with small amount of bluntness, the critical Reynolds number for the 
onset of instability increased by an order of magnitude compared to the sharp cone value. 
Gupta et al. [17] did numerical analysis for hypersonic laminar and turbulent flows 
over blunted sphere cones using time-steady viscous shock layer equations. Results were 
obtained from a spatial marching, implicit, finite difference technique, which included 
coupling of the normal momentum and continuity equations and the use of the Vigneron 
[9] pressure condition in subsonic nose region to ensure stability. 

In the present study, the PNS code of Vigneron et al. [9] was modified for the 
on-axis biconic geometry. The objective of the present investigation is to conduct a 
parametric study to determine the effect of AOA and nose bluntness on the entire flow 
field past slender bodies under different ffeestream conditions in viscous hypersonic flows. 
The slender bodies considered here are blunted cones and on-axis biconics with circular 
cross sections. Entropy was calculated at different axial locations in the streamwise 
direction for straight biconics under varying nose radii and AOA’s. By comparing 
entropy, velocity, temperature, and pressure profiles, an estimate can be made about 
entropy layer swallowing distance for the straight biconic geometry. 

The physical model and certain theoretical concepts are presented in Chap. 2. In 
Chap. 3 the theoretical formulation of the problem is provided. The solution procedure 
is explained in Chap. 4. The results obtained in the present study are discussed in detail 
in Chap. 5. Finally in Chap. 6, the conclusions are presented. 


Chapter 2 

PHYSICAL MODEL AND THEORETICAL CONCEPTS 


The calculation of hypersonic viscous flow fields past long slender axisymmetric blunt 
bodies is of prime interest to the designer of certain aerospace vehicles. A wide variety of 
flow conditions are encountered during the transatmospheric flight of these vehicles. The 
forebody of these vehicles must have finite thickness in order to sustain the pressure and 
heating rates encountered at hypersonic speeds. For selected applications, the forebody is 
modelled by spherically blunted cones and biconics. The physical model for hypersonic 
flow past a blunt slender cone is shown in Fig. 2.1. The relative positions of boundary. 


ciiuupy mid shuck layers are indicated in the figure. The first problem considered is the 

O 0 

flow past a blunt-nosed slender cone at 0 and 2 AOA. The origin of the coordinate 
system is at the virtual tip of the sharp cone with x axis along the symmetry line and y 
axis normal to it. Due* to symmetry of the flow field even at AOA, only one half of the 
flow field is computed. This case was run just to check the computer code. Numerical 
results have been obtained by employing a combination of a NS code and a PNS code. 
The results obtained are discussed in detail in Chap. 5 

O 

In the first problem, a 5 half angle cone with a nose tip radius of 0.025 m is 
selected. The second geometry considered (Fig. 2.2) is a spherically blunted, 12.84°/ 

O 0 

7 , on-axis biconic with a bend angle of 7 [10]. The maximum length of the biconic 
is 31.916 R n . The aft section starts after x=18.159 R n . For this configuration, two nose 
tips with radii of 0.00383 m and 0.04 m are considered, with the smallest nose radius tip 
approximating the sharp tip. Blunting was accomplished by keeping the fore cone angle, 
the length of the forecone, and aft section fixed, and increasing the nose radius. As the 
body dimensions are held fixed, the extent to which the bluntness effects propagate 
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downstream cannot be assessed. However, Singh and Tiwari [2, 19] have investigated 
this effect for blunted cones and ogives. Moreover, investigation in this area requires 
the minimum length of the body to be on the order of 250 R n . The bluntness effects are 
determined by comparing the solutions obtained from the blunted tip with that obtained 
from the sharp tip. 

To study the effect of AOA on the biconic configuration, the fore cone angle and 
nose radius (0.00383 m) were held fixed. Results were obtained for 8° and 12 P AOA 

O 

cases and compared with the 0 AOA case. Next, the nose radius was increased to 0.04 

O 0 

m and results were obtained for 0 and 12 AOA. 

2.1 Downstream Effects of Nose Bluntness and Angle of Attack 

In this section, the effects of AOA and nose bluntness on the flow field and wall 
quantities over the forebody of a hypersonic re-entry vehicle are considered. 

2.1.1 Aerodynamic Heating 

From the practical aspect of the design of hypersonic vehicles, the stagnation point 
heat transfer becomes a vital parameter. It has been proved that the stagnation point 
heating varies inversely with the square root of nose radius [20], i.e.. 

Hence, the leading edge heat transfer can be reduced considerably by increasing the 
nose radius or bluntness. This is the reason why the nose and leading edge regions of 
hypersonic vehicles are designed to be blunt; otherwise, the severe aerothermal conditions 
in the stagnation region would quickly melt a sharp leading edge. 

Aerodynamic heating becomes crucial at hypersonic speeds. For the flow over 
a flat plate, this is expressed as [20] 



( 2 . 2 ) 
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where C\ is the Stanton number defined by 



PooVoo{haw ^tu) 

Hence aerodynamic heating increases with cube of the velocity. By comparison, aero- 
dynamic drag is given by 

D = \ P< »VlsC d (2.3) 

which increases as the square of the velocity. Hence at hypersonic speeds, aerodynamic 
heating increases much more rapidly with velocity than drag, and this is the primary 
reason why aerodynamic heating is a dominant aspect of hypersonic vehicle design. 

A major part of aerodynamic heating for a transatmospheric vehicle is encountered 
during ascent rather than during entry; because during ascent the vehicle will accelerate 
to orbital velocity within the sensible atmosphere (using air breathing propulsion), hence 
high velocity will be combined with relatively high density, which from Eq. (2.2) combine 
to yield very high heating values. In contrast, on atmospheric entry, the transatmospheric 
vehicle will follow a gliding flight path where deceleration to lower velocities will take 
place at higher altitudes, hence resulting in lower heating rates than are encountered 
during ascent. 

2.1.2 Entropy Layer 

Nose bluntness at hypersonic speeds causes the shock wave to be detached and 
curved in the nose region. The entropy of the flow increases across the shock layer; the 
stronger the shock wave, the larger the entropy increase. Thus the streamline passing 
through the nearly normal portion of the shock will have a large entropy increase than 
a streamline passing through the weaker portion of the shock. This change in entropy 
production in the nose region generates a layer of flow with an entropy gradient, known 
as the entropy layer. Since the entropy gradients are more severe at hypersonic Mach 
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numbers, it introduces a large amount of rotationality into the inviscid flow field, as can 
be quantitatively obtained from the Crocco’s theorem 

TVs = Vh 0 - V x (V x V) (2.4) 

where Vho is the gradient of total enthalpy, and V x V is the vorticity. 

The effects of entropy gradients can be classified into two groups [21]. For two 
dimensional and axially symmetric bodies, the presence of curved shocks can produce 
large vorticity normal to the direction of velocity at the surface and, therefore, velocity 
gradients in the direction normal to the body surface. These velocity gradients change 
the boundary layer velocity profile and, consequently, the heat transfer at the surface 
of the body. In this case the boundary layer profile must be such that the velocity and 
velocity gradients outside the boundary layer must match the shear flow outside of the 
boundary layer. Both the velocity and shear flow are functions of pressure and of mass 
flow contained in the boundary layer and thus of the Reynolds number. 

The effect of entropy for three-dimensional flow can be more complicated because, 
for the three-dimensional case, the presence of entropy gradients can produce vorticity 
that has a component both normal to the stream direction and parallel to the stream 
direction of the flow outside the boundary layer. The component parallel to the stream 
direction can be large for bodies without axial symmetry and can affect heat transfer and 
boundary layer stability. 

Ferri [22] pointed out that this interaction between the rotational external flow 
and the boundary layer may in some instances invalidate the classical boundary layer 
approach, particularly when the vorticity of the external stream is of the same order as 
the average vorticity in the boundary layer. These conditions may exist, for example, 
in the combination of low Reynolds number (low boundary-layer vorticity) and high 
Mach number (high external stream vorticity because of the highly curved shock). As 
the Reynolds number is increased, the effect of external vorticity on the boundary layer 
velocity and temperature profiles is lessened. 
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The thickness of the entropy layer is a function of the leading edge bluntness, even 
a small leading edge thickness generates an entropy layer which influences the local flow 
conditions for a large distance downstream. For some distance downstream of the blunt 
nose, the boundary layer grows inside the entropy layer which eventually is swallowed 
by the boundary layer in case of the flow past blunted cones. This distance is called 
the entropy swallowing distance. For conical flows, the entropy swallowing distance is 
defined as “the location at the cone frustum where fluid, which has gone through the 
strong portion of the bow shock, has been swallowed by the boundary layer [23].” This 
is by no means a precise definition as it depends upon the shock shape, boundary layer 
assumptions, and definition of entropy layer thickness. 

The presence of entropy layer on a blunt nosed hypersonic body has an important 
effect on aerodynamic heating predictions using boundary layer techniques. In contrast, 
when the entire shock layer is treated as viscous from the body to the shock wave, 
the explicit treatment of the entropy layer is not needed. For such viscous shock layer 
treatment using a PNS approach, the interaction between the entropy layer and the shock 
layer “comes out in the wash”; no separate treatment is required, because it is contained 
within the framework of a fully viscous calculation. 

2.1.3 Bluntness and Angle of Attack Effects on Transition 

Boundary layer transition is a problem that has plagued several generations of 
aerodynamicists. The influence of nose tip bluntness on cone frustum transition remains 
an area which is poorly understood. From the results of bluntness investigations [6], it can 
be concluded that the effect of blunting on transition is very sensitive to free stream Mach 
number, with large Mach numbers producing a large rearward displacement of transition. 
The reason for this sensitivity with free stream Mach number is believed to be primarily 
related to the Reynolds number reduction associated with pressure losses across the bow 
shock. Small bluntness systematically moves the transition location rearward until the 
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maximum displacement is obtained. Large bluntness causes transition reversal. That 
is, additional increase in the nose tip radius or freestream Reynolds number produces a 
forward movement of transition. The forward movement of transition is very sensitive to 
small changes in Reynolds number or nose radius. The maximum rearward displacement 
of transition occurs under situations of low local Mach number flow. The trend of 
maximum transition displacement with free stream Mach number follows the trend of 
Reynold number reduction. Reynolds number reduction is believed to be a dominant 
effect associated with the rearward displacement of transition. 

Although transition trends on a sharp cone at AOA may defy intuition, there 
seems to be general agreement regarding the expected movement of transition. Theory 
and experiment both indicate a rearward movement of transition on the windward ray 
and a forward movement on the leeward ray. The biconic configurations investigated did 
not produce any significant delays in the boundary layer transition (as compared to cone) 
and, in fact, often promoted an earlier transition. 


Chapter 3 

THEORETICAL FORMULATIONS 
3.1 Introduction 

The design of a supersonic/hypersonic vehicle requires a detailed analysis of the 
flow field. Calculation procedures that range from simple methods to complex numerical 
techniques have been developed over years. A traditional approach to the problem is to 
decouple the flow field into an inviscid region governed by the inviscid flow relations 
and a viscous region adjacent to the surface governed by the boundary layer equations. 
Once the inviscid flow field is known, either from experimental measurements or from 
a theoretical solution, procedures can be developed to generate solutions for the thin 
boundary layer near the vehicle surface. These procedures vary in degree of sophistication 
from simple correlation to numerical programs which calculate the non-similar boundary 
layers for laminar, transitional, and/or turbulent flows. 

Two problems in the boundary layer formulation are : 

1. The uncertainties in the flow properties associated with having to determine the 
boundary layer edge required for the solution of boundary layer equations. 

2. The iterative process for locating the boundary layer edge. 

In addition, for applications where a strong interaction between the viscous and 
inviscid regions occurs, the decoupling of the flow field is no longer an acceptable 
approach. This fact is particularly important for hypersonic flow fields where the shock 
layer is relatively thin and viscous effects may influence a large portion the shock layer. 
To address these types of flow fields, procedures are developed to solve the equations 
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of motion for the entire flow field. Obviously, the NS equations may be employed to 
solve such flow fields. However, the numerical solution of the NS equations requires a 
substantial amount of computer time and storage. This is due to the many iterations which 
are typically required in solving the unsteady NS equations in time until a converged 
steady state solution is reached. Thus it is desirable to reduce the NS equations to a 
form which can be solved efficiently, while concurrently the physics of the problem is 
essentially preserved. A popular method which has proven successful for the computation 
of these classes of flow fields is the space marching of the PNS equations. In this chapter, 
the strong conservation law form of NS and PNS equations will be presented for use 
with a generalized coordinate system. 

3.2 Basic Governing Equations 

Two sets of equations are utilized to calculate the flow over blunted slender bodies. 
For the nose section, the use of the full NS equations is necessary due to the presence 
of the subsonic region; but slightly downstream of the nose region on the body, the 
flow becomes predominantly supersonic and the NS equations are simplified for space 
marching. The NS equations without body force and external heat addition can be written 
in nondimensional, strong conservation form as [24] 

dV d(E-E v ) d(F-F v ) d(G-G,) „ 

dt + Tx + Ty + ~ Wz U 1 ’ 

where 


ji 

U = {p,pu,pv,pw,pe t } 

T 

E = {pu, pu 2 + p, puv, puw, (pe t + p)u) 

rp 

E v = {0, <7 XX , Txyi T X z ) Gxx'U TxyV “H T X2 U? + (Jx } 

T 

F = {pv,pvu,pv 2 +p, pvw, {pet +p)v} 

T 

F v = { 0 , T X y Vyy, T yz , T X yU + <J yy v + Tyz w + Qy} 
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(3.2) 


P = (7 ~ l)pe 


(3.3) 


These equations are nondimensionalized with respect to freestream values [9]. The 
coefficient of viscosity p. is obtained using the Sutherland’s law 



(3.4) 
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where T re f = U S' 4 —. For perfect gas calculations, Prandtl number is taken as 0.72. 

J -L OO 

The computed entropy profiles are presented in this study. The entropy change AS is 
computed assuming a perfect gas as 


where AS = . 


AS = 


7 

7 - 1 



(3.5) 


For the spherical nose region of the body, the flow is analyzed by the code SOFIA 
[25]. The code solves the complete set of NS equations using an explicit, finite volume, 
alternating two-step scheme, and uses local time stepping to accelerate the convergence 
to steady state. It also incorporates time dependent adaptive grid to properly resolve the 
flow gradients without increasing the grid points. Once the flow field is obtained over the 
nose section, it is used to provide an initial plane solution to initialize the PNS equations. 

In the present investigation, a parabolized Navier-Stokes code developed by 


i rm 


vigneron et ai. j>j is used for analyzing the the downstream flow over the bodies. The 
PNS equations are a subset of the complete NS equations and can be used to predict many 
complex three-dimensional, steady, supersonic, viscous flow fields. These equations can 
be solved using a space marching technique instead of the time marching procedure which 
is usually employed for the NS equations. As discussed earlier, they are valid in both 
viscous and inviscid regions, and thus unlike boundary layer equations no special effort 
is needed for viscous/inviscid interactions. They can be obtained by 


1. neglecting unsteady terms. 


2. neglecting all viscous gradients in the streamwise direction. 

3. modifying the streamwise pressure gradient term in the subsonic region to permit 
stable marching of the equations from the initial data plane. 

Technically, it is not necessary to drop the unsteady terms in order to “parabolize” 
the equations (i.e., make the equations hyperbolic/parabolic in the streamwise direction). 
Some iterative space marching methods actually retain the time derivative term. 
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For space marching of these equations to be applicable, the majority of the flow field 
should be supersonic in the streamwise direction and the streamwise velocity component 
should be positive. Typically space marching of these equations will become unstable and 
fail in the presence of streamwise separation, although regions of crossflow separation 
present no special difficulties. In addition, regions of subsonic flow (immediately adjacent 
to the solid walls due to the no-slip boundary condition) should be thin, and the streamwise 
pressure gradient in these subsonic regions should not be strong. Otherwise, results could 
be poor. The following two independent variable transformations are applied to allow 
for conical effects and stretching between the body and shock 

a = x b c = — (3.6) 

x x 


£ = u »/ — if (a, 6, c) 




S’ ( L -\ 
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/a i\ 
w* ' J 


Upon combining Eqs. (3.1), (3.6) and (3.7), neglecting unsteady terms and viscous 


derivatives in the £ direction, the PNS equations are expressed as [9] 


dEi BF\ dGi 
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The flux vector Ei in Eq. (3.8) includes the pressure and, therefore, the streamwise 
pressure gradient appears in With no modification of the streamwise pressure 

gradient within the subsonic regions, the equations are elliptic in these regions, and 
thus the downstream flow field is allowed to influence the upstream flow. Therefore, a 
space marching procedure to solve the system cannot be incorporated unless this elliptic 
influence is suppressed. If a marching procedure is used, exponential growth or decay in 
the solution near the surface can occur, which in turn will cause failure of the numerical 
scheme. This failed solution is known as the departure solution. To overcome this 
difficulty, the streamwise pressure gradient is modified in the subsonic regions. One 
approach is to split the inviscid flux vector Ei as 


Ei = E\ + P 


(3.9) 


where 


rp 

E\ = {/>u, pu 2 + Clp, puv , puw, ( pet + p)u} 

P = {0,(1 -Q)p, 0,0, 0} T 

In the supersonic portion of the flow field, f!=l and no modification of the pressure 
gradient is required; but in the subsonic regions of the flow the streamwise pressure 
gradient is to be treated specially. 

3.3 Streamwise Pressure Gradient 

Numerous schemes have been introduced to handle the pressure gradient in subsonic 
regions, some of which are reviewed in this section. 

1. The obvious and simplest approximation is to drop the pressure gradient within the 
subsonic portion of the flow. This crude approximation will introduce inaccuracies in 
the solution of the flow fields where large pressure gradients are present and therefore 
has limited use. 


2. The pressure gradient is evaluated explicitly by a backward difference approximation 
evaluated at the previous plane, i.e., 


dp | _ Pi- Pi-i 

dC +1 


When this approximation is used, the stability requirement imposes a lower limit on 
the selection of the marching step size, i.e., must be larger than some (A«f) m i n 
which is provided by the stability analysis [24]. 

3. In the development of the boundary layer equations it is assumed that the normal 
gradient of pressure within the boundary layer is negligible. Consistent with this 
assumption, the streamwise pressure gradient in the PNS equations is computed at 
the first supersonic point and imposed upon the subsonic portion. This procedure 
is known as sublayer approximation. It has been observed that this approximation 
introduces instability in some cases. 

4. Earlier the flux vector Ei , which includes the pressure, was decomposed by introducing 
a parameter “ST\ If one is interested in using a space marching method for integrating 
the PNS equation set, then the inviscid eigenvalues have to be real and the viscous 
eigenvalues have to be positive and real. Based on eigenvalue stability analysis, a 
fraction Q, of the streamwise pressure gradient may be retained and evaluated implicitly 
while the remaining (1— fi) fraction is evaluated explicitly. Therefore, 


dp 


^ I implicit 4" 


(1 - explicit 


Vigneron et al. [9] determined from stability analysis that 
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where a is a safety factor, usually assigned a value of 0.8 and is the local Mach 
number in £ direction. In the regions where the flow is supersonic, > 1, ft=l, also 
< 7=1 and no approximation is incorporated. 

In the present study the fourth approximation for the streamwise pressure gradient 
has been incorporated. 

3.4 Freestream and Boundary Conditions 

The freestream and the boundary conditions are dependent on specific physical model 
considered for investigations. Discussions on typical boundary conditions are provided 
in Chap. 4, and information on specific freestream conditions is given in Chap. 5. 


Chapter 4 

METHOD OF SOLUTION 


The PNS equations are a mixed set of hyperbolic-parabolic equations in the stream- 
wise direction. Equation (3.8) is solved with a finite difference method using the Beam 
and Warming [26] algorithm. It is an implicit, non-iterative (i.e., non-iterative at a fixed 
marching station) algorithm with approximate factorization in delta form. The choice 
of an implicit algorithm is justified when the limit imposed on marching step by the 
stability condition of an explicit method is smaller than the limit required for accuracy. 
The Vigneron technique is used to suppress the numerical instabilities due to streamwise 
pressure gradients in subsonic regions. The non-iterative character of the present method 
is expected to provide better efficiency than the iterative schemes of Rubin and Lin [27] 
and Lubard and Helliwell [28]. 

For the governing Eq. (3.8) written as 


dE\ dP dF l dGi 
d( + d( + ft; + d( - 


(4.1) 


the delta form of the algorithm for constant step size A£ is [9] 
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(4.2) 


where the superscript i refers to the current known level £=i A£, i+1 is the next or 
unknown marching plane and U\ = A'Ui = U[ +1 - U[ and the derivatives dr] and 
d( are approximated with central difference operators. This algorithm has been factored 
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in terms of Ui rather than E* because the computation of the Jacobians is easier 

than the computation of as well as other reasons of computational convenience. 

The Jacobians and are evaluated in Appendix A. Since the vectors Ei, 

Fj and Gi are homogenous functions of degree one in U, the conservative form of the 
governing equations is maintained. For first order accuracy in £, the Euler implicit scheme 
is used ( 9\ = 1 and 62 — 0). The Jacobians are evaluated at level i and A e P = A t ~ 1 P. 
In the present study results were obtained with the first order scheme. If second order 
accuracy in £ is desired, one can use the Crank Nicolson scheme ( 6 \ = i, # 2 = 0) or the 
three point backward implicit scheme (9\ = 1 , 0 2 = |). 

Two approximations are made in the computation of viscous Jacobians. The 
viscosity ft is assumed independent of vector Ui and is a function of position only. The 
cross derivative terms in the ( 77 , Q plane are neglected from the Jacobians. In practice, 
Eq. (4.2) is solved in the following four steps: 

Step 1: 


'dE* 9 l A( d / dG\\ 
8U X + 1+02 dC\duJ 


A'Ui = RHS of (4.2) 


(4.3) 


Step 2: 


Step 3: 


Step 4: 
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In step 1, A'Ui represents the vector quantity 

fdEt\~ l \dEt 0iA( d (dF x \ 
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which is determined by solving the system of equations given by Eq. (4.3). This system 
of equations has the block tridiagonal structure which can be solved using the Thomas 
algorithm, where in fact, a set of uncoupled block tridiagonal systems of the type given by 
Eq. (4.3) must be solved at step 1. Once A'Ui is determined, it is multiplied by ( J 

* ( d E* \ ^ 

in step 2. As a result of this multiplication, the inverse matrix ( does not have 

to be determined in the solution process. In step 3, a second set of the block tridiagonal 

systems of equations in the tj direction are solved. Finally, in step 4, the vector of 

unknowns at station i+1 (i.e., Ui 1+1 ) is determined by simply adding A’Ui to the vector 

of unknowns at station i. The primitive variables can then be easily obtained from Ui 1+1 . 

In order to damp the spurious numerical oscillations due to central differencing in the 

cross plane, second order implicit and fourth order explicit damping terms are added. 

After adding the dissipation terms, the truncation error of the algorithm is consistent with 

the first order Euler scheme. The final algorithm is expressed (with the fourth order 

explicit smoothing [26] terms) as [9] 


dEt 
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1 + 02 V dr) _r d( J + 1 + e 2 


-A‘“ x Ei - A e P 
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j m 


(4.7) 


where V and A are the conventional backward and forward difference operators, 
respectively, and eg and ei are the coefficients of explicit and implicit dissipation. 

Due to symmetry of the flow field, only one half of the flow field is computed. 
The generalized coordinates rj and £ are defined such that the computational plane has 
a square shape of side unity with uniform spacing in both directions. Therefore, the 
correspondence between physical and computational plane is, for 1 < j < NJ and 1 < 
k < NK, as follows 
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where 


T) = (k — 1 )Atj 

C = (i - 1)AC 
y{h k) = VBii) + 

z(j \k) = ZB(j) + s(hj,P) sin [<f>(j)} 


At/ = 


1 

NK - r 


AC = 


1 

NJ - 1 


The variable s is the stretching function depending upon rj, 6 and a stretching parameter 
/S, and is expressed as 


sihjiP) = 



(4.8) 
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grid points are clustered near the body to properly resolve the gradients. In the present 
study, a grid of 30x50 points was used, with 30 points in the circumferential direction, 
and 50 points in the normal direction. Figure 4.1 shows the grid for the biconic geometry 

o o 

at 0 AOA and Fig. 4.2 shows the grid at 12 AOA for the same geometry. 


4.1 Boundary Conditions 

Boundary conditions are to be specified at the body surface and at the outer boundary 
which in the present study is a bow shock. At the body no slip, zero pressure gradient 
and constant wall temperature conditions are imposed. These equations are analyzed 
mathematically in this section, following the procedure outlined in [29]. 

One may note that 

(Ape t ) u = (/>e<) i+1)1 - (pe t ) itl (4.9a) 

At the surface, pet = pe because the velocity at the wall is zero. In addition. 
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Fig. 4.2 Grid for 12.84°/ 7° biconic configuration at 12° AOA. 
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For a uniform surface temperature, T w (, +1 ) = = r w . Therefore, 

(A/>e*)i,l = - Pi, l) 

where <}>' is defined as <j>' — ^zfysp • Equation (4.10a) may be rearranged as 
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By definition, one may also write 


( A P)i,i = P(i+ 1,1) “ P(i, 1) 
Equation (4.10c) is subtracted from Eq. (4.11) to yield 


(4.11) 


( A P )i,l “ ti^P^i,! ~ P(»'+ 1,1) “ P(i,l) ~ (p(i+ 1,1) - P(i,l)) ~ 0 


. Thus, the required boundary conditions at the body for constant surface temperature 


can be summarized as 


— ^( A Pet)i t i — 0 


(Apu),-i = 0 
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(A / 9u)- 1 = 0 (4.12c) 

(A/ntj), j = 0 (4.1 2d) 


(Apet), ! - V»i+ 1 r (ti< >2 + w? 2 + w, ? , 2 )(Ap) f>2 - « it2 (Apu) i>2 


t>i, 2 (A/w) ,- 2 - u>,-, 2 (Apu>) <>2 + (Ape<) - 2 


(4.13) 


where * +1 = = (£). +i - 


At the outer boundary, a shock fitting technique is used and the flow variables 
behind the shock are calculated from the Ranldne-Hugoniot relations. The plane of 
symmetry boundary conditions are computed by reflection and they are imposed explicitly 
as well as implicitly. In addition to the boundary conditions, initial conditions are also 
necessary. As mentioned earlier, near the nose tip the initial plane solution to initiate the 
marching procedure is provided by the Navier-Stokes equations. 


4.2 Shock Fitting Procedure 

The bow shock generated by a slender body in a supersonic/hypersonic flow field may 
be selected as the outer boundary of the domain and determined as a part of the overall 
solution. Obviously this procedure has some advantages compared to an outer boundary 
set at the free stream where the bow shock must be captured. First, the number of grid 
points in the domain may be decreased. This is due to the fact that additional grid points 
must be located in the free stream if the bow shock is to be captured. Second, shock 
smearing, which is a consequence of shock capturing, will not appear in the solution 
when shock fitting is used. Of course, these advantages are accompanied by some 
disadvantages. Most shock fitting procedures are explicit, and therefore some additional 
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stability requirements are imposed. Furthermore, additional sets of equations must be 
used to determine the shock location. Usually the Rankine-Hugoniot relations are used 
for this purpose. 

The initial conditions at “i=l” provide all the required data including the shock 
slope. To generate the grid at the next station, i.e., at “i=2, the shock is linearly 
extrapolated. Therefore, the grid system at “i=2” is generated and the PNS equations are 
solved for j=l (at the surface) to j=jmax (at the shock, i.e., jmax is the grid point just 
behind the shock). The finite difference form of the PNS equations are modified and a 
one sided difference approximation is employed at the shock. Once the PNS equations at 
“i=2” are solved, all the flow properties including the pressure are known. However, the 
shock location at “i=2” was extrapolated from the previous station. Therefore an updating 
procedure must be used to modify the shock properties and compute a new shock slope 
at “i=2”. The procedure is described in two steps. First, an equation for the shock shape 
is derived, and subsequently the Rankine-Hugoniot equations are introduced. 

The shock standoff distance 6 is obtained from the values at £ through an Euler 
integration as [9] 


t(( + Aft = f({) + -A{ (4.14) 

The problem is to determine the slope || at the station £. The inward unit normal to 
the shock is given by 


N,= 
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and the derivatives with respect to £ are taken along the shock. If Vi denotes the 
upstream flow velocity normal to the shock, then 

2 


Vi 2 = (jVVoo)' 

Upon substituting for Voo and N s , this equation can be solved for §| as 
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Wl + Vlyjtih + tl 


(t 2k) 


dn 

di 


(4.16) 


(4.17) 


where 


. ( dz 

h ~ v<x V dC I ~ W °° 

*>-*.-* 


( 

\d() 


h - (I) + (I) 

i4= c o^(|)- S in^(||) 

Uoo = cos a, Voo — ~ sin a, tOoo = 0 

The normal component of velocity behind the shock, Vi and the density, p\ are 
expressed by (Rankine-Hugoniot relations) 


Vi = 


7+1 fp + 7-1 




7 + 1. 


(4.18) 


Pi = 


i+p {&) 


(4.19) 


where. 


P = 


The pressure behind the shock Pi can be calculated by the application of a one-sided 
version of the finite difference algorithm, once the new shock position is determined. 
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The metric coefficient i s obtained by differentiating Eq. (4.8 ) 


dn _ 2/3 s_ fdS _ S\ 

9a {^-( £ f 2 ) 2 }'"(f±r) p ^ 


(4.20) 



Chapter 5 

RESULTS AND DISCUSSION 

The numerical schemes described in Chap. 4 were used to obtain results for two 
different geometries, slender blunted cones and on-axis biconics. All calculations were 
performed under the assumption of a laminar flow and perfect gas. For the forebody, the 
flow is calculated by the parabolized Navier-Stokes equations with the starting solution 
for the nose section provided by the full Navier-Stokes equations. The parametric studies 
were conducted to determine the effect of nose bluntness and AOA on surface quantities 
and flow field variables for straight biconics. The effect of AOA on surface quantities 
and shock stand-off distance for a blunted cone was also investigated to check the PNS 
code. The freestream conditions for the blunted cone and on-axis biconic considered for 
this study are given in Tables 5.1 and 5.2 respectively. These conditions were obtained 
from [10]. 

The results for slender blunted cone are presented first and this is followed by 
the results for the straight biconic. The results obtained for the biconic configuration are 
subdivided into several subsections, and these are presented in a somewhat logical order. 

5.1 Results for Slender Cone 

The effect of nose bluntness on flow field and surface quantities was investigated 
extensively in [2, 19] for blunted slender cones and ogives at a=0°. Here , selected 
results obtained for a=0° and 2° are presented to demonstrate the influence of AOA on 
different quantities. 

Figures 5. 1-5.4 show the variation of surface quantities with the axial distance 
for a 5 cone at 0 and 2 AOA for freestream Mach number 20. The Reynolds number 
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Table 5.1 Freestream flow conditions for 5° cone, a= 0° and 2° 
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Quantity 

Me =20 

Poo , N/m 2 

171.0 

Twall» k 

1000 

Too,k 

261.3 

Pr 

0.72 

Re(R n ) 

22303 

R„, m 

0.025 

7 

1.4 


Table 5.2 Freestream flow conditions for 12.87°/ 7° 


° ° 0 
bicone; a = 0 , 8 and 12 


Quantity Moo = 6.89 


'oo >N/m 2 

2182 

T W all» k 

300 

To© » k 

1604 

Pr 

0.72 

Reoo 

4.5x1 0 5 

R n , m 

0.00383 

7 

1.3 



36 


Re (R n ) considered is based upon the nose radius. The nose radius considered is R n =0.025 

O O 

m. A ray at <f>=Q is considered to be most windward and the one at </>=180 is considered 
to be most leeward. 

The surface pressure distributions are presented in Fig. 5.1. The pressure has been 
non-dimensionalized with respect to the freestream value. The x distance is measured 
from the virtual tip of the cone. The results are plotted from the tangency point onward. 

O 

Consider the incoming flow at 0 incidence, in the stagnation region, the pressure is 
very high but as the flow moves downstream, the shock strength decreases, thereby 
decreasing the pressure. The high pressure flow continues to expand along the cone 
frustum producing a favorable pressure gradient as evident from the figure. The favorable 
pressure gradient region for present conditions extends to about x=60 R n . However, by 
now the flow has overexpanded and must recompress in order to attain the conical wall 
pressure value far downstream of the nose. This recompression then results in an adverse 
pressure gradient region. The extent of overexpansion and recompression is a function 
of bluntness and freestream conditions. The surface pressure reaches the conical value 
at about x=160 R„. 

O 

At Mach 20, even if the incoming flow is at small incidence a=2 , it has significant 
effects on surface pressure on windward and leeward planes. For the windward plane the 
favorable pressure gradient extends to about x=48 R n . The pressure attains a constant 
value at about x=188 Rn> which is approximately 66% higher than the conical value 
at a=0°. On the leeward plane, the favorable pressure gradient extends to about x=84 
R n . The overexpanded flow compresses very slightly and almost becomes constant at 
approximately x=136 R„. This is about 63% lower than the conical value at a=0°. On the 

O 0 

sideward meridian <f >= 90 , the surface pressure behaves in the same fashion as at 0 AOA. 

Figure 5.2 shows the variation of shock stand-off distance as a function of the 
axial distance. The shock stand-off distances are measured from the body surface. The 

O 

shock stand-off distance is nondimensionalized with respect to nose radius. At 0 AOA, 



o=2, <£=180° 

a=0* 

a=2\ 0=0° 

a=2’,<f>=90° 



7 = 1.4 
= 20 

Rn=0.025 m 



X, m 

Fig. 5.2 Variation of shock standoff distance with 

• 0 0 
axial distance for blunted cone at 0 and 2 AOA. 




Fig. 5.3 Variation of wall heat transfer with axial 
distance for blunted cone at 0° and 2° AOA. 






Fig. 5.4 Variation of skin friction coefficient with 
axial distance for blunted cone at 0° and 2° AOA. 
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the point of inflection occurs at about 170 nose radii downstream. At 2° AO A the shock 

gets displaced towards the leeward meridian. On the leeward side the shock is totally 

curved and inflection point disappears. On the windward side the shock gets compressed 

and the measured stand-off distances fall below the 0 AOA case. On the windward 

meridian the inflection point is about 96 R n downstream. On the sideward plane (<p=90 ) 

0 

the shock shape lies very close to 0 AOA curve. 

The variation of wall heat flux with axial distance is illustrated in Fig. 5.3. The 
wall heat flux has been nondimensionalized by pooV^o- The heating rates are very high 
in the nose region and decrease continuously over the afterbody due to lower pressure. 
Because of the favorable and adverse pressure gradients near the nose region, the heat 

flux variation show a similar behavior to that of the surface pressure. Even at small 

0 

incidence (a=2 ), the wall heat flux increases on the windward plane considerably. On 
the leeward side there is a large decrease in the heat flux in comparison to 0° AOA case. 

O 0 

On the sideward plane (</>=90 and a= 2 ) the heating rates are little higher as compared 

O 

to the a=0 case. 

O 0 

Variation of skin friction coefficient with axial distance at 0 and 2 AOA is 
shown in Fig. 5.4. It follows the same general trend as the wall heat flux. This similarity 
can be explained qualitatively on the basis of Reynolds analogy. The skin friction on 
the leeward side decreases by an order of magnitude compared to the windward side at 

O 

a=2 . The reason is that , on the leeward side, the entropy and boundary layers are very 
thick, and gradients of velocity near the wall are small. 

5.2 Results for the Straight Biconic 

The various results obtained for 12.84°/ 7°, straight biconic configuration are pre- 
sented in this section. All the results were obtained for a freestream Mach number of 
6.89, 7=1.3, cold wall condition (^- = 0.187) and the free stream Reynolds number 
per unit length, Re,*, = 4.5 x 10 5 (Table 5.2). For numerical calculations, the Reynolds 
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number based on nose radius was considered. Two nose radii were considered to study 
the bluntness effects; R n =0.00383 m and 0.04 m. The AOA was varied from 0° to 8° 
and 12° to study its influence on the flow field and surface quantities. 

5.2.1 Experimental Verification 

A few important results obtained by PNS code are validated by the experimental 
results of Miller et al. [10] for the same geometry and freestream conditions. Figure 5.5 
shows the variation of wall heat flux Q w (MW/m 2 ) as a function of nondimensionalized 
axial distance at a=0°. The experimental result is accurately predicted by the PNS code. 
Figure 5.6 shows the variation of Q w with x/L at 12° AOA. The PNS code slightly 
underpredicts the windward and leeward heating. This is due to the fact that at a=12°, 
the boundary layer is thinner on windward side which requires more number of grid 
points to accurately model the temperature distribution through the boundary layer. Also 
the flow separates on the leeward side. 

5.2.2 Nose Bluntness Effects 

The effect of nose bluntness on wall quantities at 0° and 12° AOA are presented 
in Figs. 5.7-5.14. The results shown in Figs. 5.7-5.11 correspond to a=0°. Figure 5.7 
shows the variation of nondimensionalized surface pressure, ■£“- with dimensional axial 

* oo 

distance, X in meters. There is no appreciable change in surface pressure even if the 
nose bluntness is increased by one order of magnitude for the present biconic geometry. 

In Fig. 5.8 the effect of nose bluntness on wall heat flux Q w (MW/m 2 ) is shown. 
If the nose bluntness is increased by an order of magnitude, heating on the fore-cone 
section decreases four times and approximately two times on the aft-cone section. This 
behavior was also observed for the skin friction coefficient which is shown in Fig. 5.9. 

In Fig. 5.10 the variation of dimensional shock stand-off distance 6 , in meters as 
a function of dimensional axial distance X, in meters at R n =0.00383 m and 0.04 m are 
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shown. Blunting of the sharp nose, causes an increase in the shock stand-off distance 
by an order of magnitude. 

Figures 5.11-5.14 correspond to q= 12°. The surface pressure variation for 
leeward and windward meridian is shown in Fig. 5.11 and the result looks similar 
to a=0° case. The bluntness effect on fore-cone and aft-cone heating is presented in Fig. 
5.12. The fore-cone heating is more sensitive to nose bluntness than aft-cone heating. 
Heating on the most windward ray decreases by four times for the fore-cone section and 
about two times for the aft-cone section. Heating on the most leeward ray decreases 
by five times for the fore-cone section, but there is a little change for aft-cone heating. 
The bluntness effect on skin friction coefficient shown in Fig. 5.13 is analogous to wall 
heat flux. In Fig. 5.14, the dimensional shock shape r S h, m is plotted as a function of 
dimensional axial distance X, m. The two sets of curves correspond to R n =0.00383 m 
and 0.04 m, respectively. The distance. r s}j is measured from the axis of the biconic. 
The leeward and windward shock shapes increase by an order of magnitude with the 
increase in nose bluntness. 

O 

Figures 5.15-5.17 represent the circumferential variation of wall quantities at 0 , 

O 0 

8 , and 12 AOA. These curves are plotted for nose bluntness, R n =0.00383 m and axial 
location x/R n =18.016. Figure 5.15 shows the variation of the coefficient of pressure, Cp, 
with meridian angle <j>. It is seen that the value of C p increases on the windward side and 

decreases on the leeward side with increase in AOA. The value of Cp is maximum on the 

0 0 

windward side (<£=0 ) and decreases monotonically towards the leeward side (<£=180 ). 

O 

At 0 AOA, C p remains constant in the circumferential plane. The circumferential plot 
of wall heat transfer (Fig. 5.16) and skin friction coefficient (Fig. 5.17) exhibit a similar 
behavior as the coefficient of pressure. 

In Figs. 5.18 and 5.19, the circumferential variation of Cp and wall heat transfer 
are shown at 0°and 12° AOA for R n =0.04 m and axial distance x/R n =30.89. In Fig. 5.18 

O O 

the value of Cp starts increasing slightly on the leeward side between <£=150 to <£=180 . 




X 

L 


Fig. 5.5 Comparison of experimental and numerical results of variation 
of wall heat flux with axial distance for straight biconic at 0° AOA. 
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Fig. 5.6 Comparison of experimental and numerical results of variation 
of wall heat flux with axial distance for straight biconic at 12° AOA. 
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Fig. 5.9 Effect of nose bluntness on skin 
friction coefficient for on-axis biconic at 0° AOA. 


| 









1.6E-002 



Fig. 5.13 Effect of nose bluntness on skin friction 

O 

coefficient for on-axis biconic at 12 AOA. 
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7= 13 
Moo = 6.89 



Fig. 5.16 Circumferential variation of wall heat transfer 
for on-axis biconic at 0°, 8° and 12° AOA, x=18.016 R n . 




Fig. 5.17 Circumferential variation of skin friction coefficient 
for on-axis biconic at 0°, 8° and 12° AOA, x=18.016 R n . 
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Fig. 5.19 Circumferential variation of wall heat transfer for 
on-axis biconic at 0 and 12 AOA, R n =0.04 m, x=30.89 R, 
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The circumferential plot of Q w (Fig. 5.19), follows the usual monotonic trend upto 

O O 

</>=140 , after that there is a sudden decrease in the value, and between </>=150 to 

O 

</>=180 the value of Q w starts increasing. As mentioned earlier this peculiar behavior is 
observed due to flow separation and formation of vortices on the leeward side. The flow 
separation in the cross flow plane can be observed clearly in Fig. 5.20 which shows the 
velocity vector plot in the circumferential plane. 

5.2.3 Effects of Angle of Attack 

Figures 5.21-5.30 show the variation of surface quantities with axial distance for on- 
0 0 

axis biconic (12.84 / 7 ) at different angles of incidence for free stream Mach number 
6.89. In Figs. 5.21-5.25 the nose radius is 0.00383 m. Figure 5.21 shows the surface 

0 0 o 

pressure distribution at 0 , 8 and 12 AOA. The surface pressure is nondimensionalized 

by the free stream static pressure and plotted as a function of x/L. The variation of the 

windward surface pressure with x/L for the on-axis biconic illustrates the overexpansion 

of the flow from the spherical nose to the fore-cone section and the expansion of flow 

from fore-cone to aft-cone section. The region of this overexpansion-recompression of 

the windward flow from the nose diminishes with increasing AOA (effective cone half 

angle, which is equal to 0 C + a). The windward surface pressure distribution increases 
0 

with AOA. At a=8 the surface pressure is nearly constant for 0.36 < x/L < 0.57 and 

O 

at a=12 for 0.28 < x/L < 0.57. At the fore-cone aft-cone junction the pressure drops 
suddenly due to discontinuity, expands a little and then attains almost a constant value. 

O O 

At 8 and 12 AOA the surface pressure on the leeward side decreases monotonically on 
the fore-cone and aft-cone regions. With increasing AOA the surface pressure distribution 
on the leeward side decreases as is evident from the figure. 

Figure 5.22 shows the variation of wall heat transfer (MW/m 2 ) distribution with 
the axial distance which has been nondimensionalized by the total length of the biconic at 
0°, 8° and 12° AOA. Heating on the windward ray increases with AOA for the straight 
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biconic geometry. Leeward heating for the on-axis biconic decreases as a is increased. 

O 

However, at a=12 on the aft section the heating rate first decreases and then slightly 
increases. This can be attributed to the flow separation at high AOA and formation of 
vortices. 

The variation of skin friction coefficient with nondimensionalized axial distance 

. 0 0 o 

is shown in Fig. 5.23 for the straight biconic at 0 , 8 and 12 AOA. Skin friction 
coefficient behaves in the same way as wall heat transfer. 

Figure 5.24 shows the variation of shock stand off distance with nondimensional- 
ized axial distance x/L for the on-axis biconic at various angles of attack. As it is evident 
from the figure that, an inflection in the shock downstream of the spherical nose occurs 

O O 

on the windward side at 8 and 12 AOA. With increase in AOA the point of inflection 
shifts towards the nose section. The shock standoff distance on the leeward side is more 
sensitive to AOA than on the windward side. For straight and bent biconic geometry. 
Miller et al. [1] has given an empirical relation to measure shock standoff detachment 
distance on leeward side at the junction of aft-fore cone sections. On the leeward side 

O 

the inflection point disappears. At 0 incidence the shock standoff distance increases 
more rapidly after the fore-cone aft-cone junction. This behavior is also observed on the 
windward side at 8° and 12° AOA. However, at 12° AOA, the shock standoff distance 
reach almost a constant value far downstream of the aft section. On the leeward side 
the shock stand off distance increases more linearly as the flow is marched downstream 
of the nose region. 

O O 

Figure 5.25 shows the shock shape at 8 and 12 AOA for the straight biconic 
geometry. From the figure it is evident that on the leeward meridian the shock shape is 
more sensitive to AOA as compared to the windward meridian. 

Figures 5.26-5.30 show the variation of surface quantities for the on-axis biconic 

O O 

at 0 and 12 AOA for R„ =0.04 m. In Fig. 5.26 the surface pressure distribution is 
shown, which looks similar to Fig. 5.21. In Fig. 5.27 the variation of wall heat flux 
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Fig. 5.22 Variation of wall heat transfer with axial 
distance for on-axis biconic at 0°, 8° and 12° AOA. 
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(MW/m 2 ) has been plotted as a function of nondimensionalized axial distance x/L. At 

0 

12 AOA, the leeward heating shows peculiar behavior. In the aft-section of the biconic, 
the heat flux increases as we march downstream. As mentioned earlier, at high AOA 
the leeward flow is characterized by flow separation and formation of two symmetrical 
counter rotating vortices that reattach on the most leeward ray [10]. At low Reynolds 
number, the leeward flow should be free of any secondary vortices [30]. This vortex 
reattachment results in the higher heating on the most leeward rays. Heating rates start 
increasing slightly on the windward side far downstream of the aft-cone section. The 
variation of skin friction coefficient with axial distance is presented in Fig. 5.28. It 
follows the same general trend as does the wall heat transfer. 

Figures 5.29 and 5.30 show the variation of shock stand-off distance and shock 
shape as a function of axial distance x/L for the on axis biconic. The results are similar 
to one shown in Figs. 5.24 and 5.25, respectively, but with the increase in hluntness, the 
shock stand-off distance and shock shape increase by an order of magnitude. 

5.2.4 Effects of Nose Bluntness and Angle of Attack on Flow field Quantities 

Figures 5.31-5.33 show the variation of streamwise velocity as a function of 
normal distance r for Mach 6.89 perfect gas flow over straight biconic at different axial 
distances. The velocity has been nondimensionalized with respect to its free stream value. 
It should be noted that since the shock fitting procedure is used, the normal profiles are 
plotted within the shock layer only. There is a significant difference between the profile 
shape near the nose and that which is seen far downstream. The thickness of the boundary 
layer increases as the flow is marched in the axial direction downstream. In Fig. 5.31 the 

O 

streamwise velocity profiles are shown at 0 AOA, where as Fig. 5.32 shows the velocity 

O 

profile at 12 AOA on leeward side. By comparison we find that in Fig. 5.31 the slope 
of the profile changes on the outer edge of the boundary layer between x/R n =7.04 and 
x/R„ =18.73. This causes an appreciable influence on the wall quantities. Also, in Fig. 
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5.32, the presence of large velocity gradients at distances far downstream indicates that 
boundary layer is relatively thicker as compared to Fig. 5.31. Figure 5.33 shows the 
velocity profiles on windward side for 12° AO A. At any particular axial location the 
boundary layer is thin on the windward plane compared to the leeward plane. The shape 
of the velocity profiles also looks different in comparison to the leeward plane. This is 
because the corresponding entropy layer is very thin. 

O 

Figure 5.34 shows the temperature profiles for the straight biconic at 0 AOA. The 
temperature has been nondimensionalized with respect to its free stream value. These 
profiles have been plotted for the cold wall condition. The temperature starts rising 
as we move from the wall towards the bow shock because of the frictional heating or 
strong viscous dissipation within the thermal boundary layer. The peak temperature 
occurs in the thermal boundary layer. Although there is a temperature jump behind 
the bow shock, this increase in temperature is not so intense compared to the peak 
temperature in the boundary layer. Therefore, just outside the boundary layer or entropy 
layer the temperature starts decreasing until the bow shock is approached. The bow 
shock is stronger near the nose region and its strength decreases farther downstream. 
Correspondingly, the peak temperature observed is higher near the nose region compared 
to that at axial distances farther downstream. However, on the aft-section of the biconic, 
the shape of the temperature profiles at the outer edge change significantly due to the 
presence of a thick entropy layer. Figures 5.35 and 5.36 show the temperature profiles at 

O 

12 AOA for the leeward and windward meridians respectively. On the leeward plane, 
strong temperature gradients are observed at axial distances far downstream from the nose 
region. The temperature profiles become almost parallel to each other, a short distance 
downstream of the nose section. The thermal boundary layer thickness in the leeward 
plane is almost twice that which observed in the windward plane on the aft-section of the 
biconic. Because of a thin boundary layer in the windward plane, the peak temperature 
close to the nose region is high in comparison to what is observed in the leeward plane. 
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Figure 5.37 show the entropy profiles for the straight biconic at 0 AOA. Entropy 
change was computed according to Eq. (3.5). Figures 5.38 and 5.39 show the entropy 

O 

profiles at 12 AOA for the leeward and windward planes respectively. These profiles 
follow the same general trend as the temperature profiles. Variation of pressure in the 
normal direction is experienced mostly in the vicinity of the nose region, and very close to 
the wall. As the flow is marched downstream in the axial direction, this effect diminishes. 
Hence in computing entropy change the flow field temperature has a dominating effect 
compared to the pressure. By comparing entropy and velocity profiles we can conclude 
that the entropy layer is not swallowed until x=30 R n . For calculating the swallowing 
distance it is necessary to march further downstream, which requires the body geometry 

0 O 

to be modified. Pressure profiles for the straight biconic at 0 and 12 AOA are shown in 
Figs. 5.40-5.42 respectively. The pressure has been nondimensionalized with respect to 
the freestream Dressure. Onlv in the nnse region are the nressnre crrariipnts strnncr ( close 

i * o ' — * r o ■ _ ~ " o ' 

to the wall), but as the flow is marched progressively towards the fore-cone section, the 
magnitude of these gradients decreases and Pitot pressure almost becomes constant in the 
normal direction. The magnitude of pressure on the leeward side is much less than the 
windward side, as can be seen clearly from Figs. 5.41 and 5.42. 

0 

Figures 5.43-5.46 show the contour plots for the on-axis biconic at 0 AOA. 
These contours are plotted within the body surface and bow shock at streamwise axial 
distance x=11.9 R n and for nose radius R n =0.00383 m. Mach contours are shown in Fig. 
5.43. The Mach number continuously decreases from the bow shock towards the body. 
Figure 5.44 shows the density contours. There is very slight change in the density within 
the shock layer. The reason is clear from the pressure contours shown in Fig. 5.45. 
Temperature contours are plotted in Fig. 5.46. Their behavior is clear from temperature 
profiles shown in Fig. 5.34. 

Figures 5.47-5.50 show the contour plots at 12° AOA. These are plotted at axial 
location x=26 R n for R n =0.00383 m. Figure 5.47 shows the Mach contours. From 
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the figure it is clear that the shock is displaced on the leeward side and is no more 
concentric with the body. A viscous ’hump’ appears on the leeward side which indicates 
that boundary layer is quite thick [3]. 

On the windward side the value of Mach number is low compared to that on the 

leeward side because the corresponding boundary layer is thin. In hypersonic regime, the 

compressible boundary layer thickness is directly proportional to the square of the Mach 

number. In the thermal boundary layer (Figs. 5.35 and 5.36) the temperature is very 

high. This affects the local velocity of sound and hence the Mach number. The density, 

pressure, and temperature contours are plotted in Figs. 5.48-5.50, respectively. Figures 

5.51-5.54 show the contour plots of Mach number, density, pressure, and temperature 
0 

respectively at 12 AO A and R n =0.04 m. These contours are also plotted at x=26 
R n . They look similar, as observed for the case of R n =0.00383 m but physically they 
correspond to a different axial location, since the hlnntness has been increased by one 
order of magnitude. 
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Fig. 5.32 u velocity profile at leeward side for on-axis biconic at 12 AOA. 
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Fig. 5.33 u velocity profile at windward side for on-axis biconic at 12° AOA. 












Fig. 5.35 Temperature profiles at leeward side for on-axis biconic at 12° AOA. 












Fig. 5.38 Entropy profiles at leeward side for on-axis biconic at 12° AOA. 
















Fig. 5.42 Pressure profiles at windward side for on-axis biconic at 12° AOA. 


























Chapter 6 
CONCLUSIONS 


The PNS code is successfully applied to study the effects of leading edge bluntness 
and angle of attack on the flow field and wall quantities of a spherically blunted on-axis 
biconic configuration in air at Mach 6.89 and 7=1.3. The variation of wall quantities at 

O O 

low incidence (a=2 ) and high Mach number (M**, = 20) for a 5 spherically blunted 
cone is also studied. The downstream influence of nose bluntness for the blunted cone 
persists up to 200 nose radii approximately. The complete flow field is accurately 
predicted using a combination of the full Navier Stokes equations in the nose region 

O O 

and parabolized Navier Stokes equations farther downstream. For the on-axis 12.84 / 7 
biconic configuration, the influence of entropy layer generated by the shock curvature on 
the flow field is also addressed. It is found that the entropy layer affects the development 
of the boundary layer far downstream of the nose section. The entropy layer is not 
swallowed by the boundary layer for the straight biconic geometry because the flow 

solution is marched only 32 nose radii downstream of the nose region. 

0 0 

With the increase in AO A from 0 to 12 , heating on the windward ray increases 
for the blunted on-axis biconic, and decreases on the leeward ray. The same trend is 
observed for the surface pressure and skin friction coefficient, but shock stand-off distance 
increases on the leeward side and decreases on the windward side. A similar reasoning 
can be applied for the results obtained for the spherically blunted cone. The decrease 
in the leeward side heating rate is more sensitive to AOA compared to the increase in 
the windward side heating. Heating on the windward side of the on-axis aft-cone is 
roughly half of that observed on the fore-cone, and leeward side heating is one order of 
magnitude less than that on the windward side. For the on-axis biconic, the heating on 
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the leeward side at the aft-cone section starts increasing a little at a=12°. This effect is 
more pronounced when the nose bluntness is increased by an order of magnitude. This 
trend is attributed to flow separation on leeward side and formation of vortices. The 
circumferential plot of wall heat transfer and coefficient of pressure and also the cross 
flow velocity vector plot show this effect more clearly. 

With the increase in nose bluntness by an order of magnitude, the windward 
heating for the on-axis fore-cone section decreases by a factor of four, and by a factor 

O 

of two on the aft-cone section at 12 AOA. Heating on the leeward ray also decreases 
by a factor of about five for the fore-cone section, but there is a little change for aft- 
cone section. A similar trend is observed for the skin friction coefficient. The surface 
pressure distribution is not sensitive to leading edge bluntness. The shock shape for the 
blunted bicone looks similar to the sharp one, but the shock radii increase by an order 
of magnitude. 
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APPENDIX 


CALCULATION OF JACOBIANS 

The Jacobians and are given by 
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These Jacobians have an inviscid part and a viscous part 
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The inviscid part can be written as a linear combination of 
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The viscous part of the Jacobians is 
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where 
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and (.)rj indicates derivative with respect to 77. Similarly 
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where 
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ani ^ (.)r indicates derivatives with respect to In these viscous Jacobians, the cross 
derivative viscous terms have been neglected and the coefficient of molecular viscosity 
has been assumed to depend only on the position, not on the vector U. Finally, the 
Jacobian is given by 
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